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Abstract 

Recent progress in the development of efficient computational algorithms to price 
financial derivatives is summarized. A first algorithm is based on a path integral 
approach to option pricing, while a second algorithm makes use of a neural network 
parameterization of option prices. The accuracy of the two methods is established 
from comparisons with the results of the standard procedures used in quantitative 
finance. 
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1 Introduction 

The development of efficient computational algorithms to price and hedge 
financial derivatives is a particularly lively topic of the research work in quan- 
titative finance and econophysics. 

In the classical theory of option pricing [1] by Black and Scholes, and Merton 
(BSM), the price of a financial derivative is given by a deterministic partial dif- 
ferential equation. In principle, by solving this equation with the appropriate 
boundary conditions, the value of the derivative of interest can be obtained. In 
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practice, the BSM equation can be analytically solved only in the most simple 
case of a so-called European option. Actually, if we consider other financial 
derivatives, which are commonly traded in real markets and allow anticipated 
exercise or depend on the history of the underlying asset, analytical solutions 
do not exist and numerical techniques are necessary. As discussed in the liter- 
ature [2], the standard numerical procedures are binomial trees, the method of 
finite differences and the Monte Carlo simulation of random walks. Since bino- 
mial trees and finite difference methods are difficult to apply when a detailed 
control of the paths is required, the Monte Carlo simulation is the mostly used 
method to price path-dependent options. However, Monte Carlo is known to 
be time consuming if precise predictions are required and appropriate variance 
reduction techniques must be introduced to save CPU time [2]. 

This paper summarizes recent progress in the development of novel computa- 
tional methods to price options, as alternatives to traditional numerical pro- 
cedures used in finance. A first algorithm is based on a path integral approach 
to option pricing, as described in detail in [3] . A second algorithm relies upon 
neural networks and represents work in progress [4]. 



2 The path integral algorithm 

The path integral method is an integral formulation of the dynamics of a 
stochastic process [5,6]. It allows to compute the transition probability associ- 
ated to a given stochastic process. We consider, for definiteness, as stochastic 
model for the time evolution of the asset price S the standard BSM Brownian 
motion, driven by the stochastic differential equation [1,2] d(lnS) = Adt+adw, 
where a is the volatility, A = (fi — ^ 2 /2), ji is the drift parameter and w is a 
Wiener process. The probability associated to the above process for a transi- 
tion from Zi = In Si to Zf = In 5/ over a time interval At is given by 



If we consider a finite-time interval [t', t"\ and we apply a time slicing in terms 
of n + 1 subintervals of length At = (t" — t')/n + 1, the finite-time transition 
probability can be expressed as a convolution of the form [3] 




(1) 



+oo +oo 



1 



1 



n+1 



*dz\ • • • dz 1 



n 



^(27ra 2 At) n + 1 



cxp 



2a 2 At 



k=l 



X)[z fc -(z fc _i + AAt)n.(2) 



— oo — oo 



2 



Once the transition probability is known, the price of an option can be com- 
puted as the conditional expectation value of a given functional of the stochas- 
tic process, involving integrals of the form [3,6] 

EPilSi-i] = J dzip^z^O^). (3) 



For example, for an European call option at the maturity T the quantity of 
interest will be max {St — X, 0}, X being the strike price. 



2.1 Transition probability 



By means of appropriate manipulations of the integrand entering eq. (2), it 
is possible to obtain an expression for the transition probability suitable for 
an efficient numerical implementation [3]. If we define y^ = z k — kAAt, k = 
1, ... ,n, we can rewrite eq. (2) as 
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in order to get rid of the drift parameter. A convenient quadratic form can 
be extracted from the argument of the exponential function, to arrive, after 
some algebra, at the following formula for the probability distribution [3] 
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where we have introduced new variables hi obeying the relation 
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The gaussian variables hi can be extracted from a normal distribution with 
mean (yoOu + y n +iO n i) 2 /mi and variance a 2 At/mi, where TOj are the eigen- 
values of a real, symmetric and tridiagonal matrix M and are the matrix 
elements of the orthogonal matrix O which diagonalizes M, with Wi = O^yj. 
The probability distribution, as given by eq. (5), is an integral whose kernel 
is a constant function and can be very efficiently and precisely computed as 
shown in [3]. 
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2.2 Option pricing 



In the case of options with possibility of anticipated exercise before the expira- 
tion date, it is necessary to check at any time ti, and for any value of the stock 
price Si, whether early exercise is more convenient with respect to holding the 
option for a future time. Therefore, the value of the option at a time slice i is 
given by 

Oi(Si) = max {Of (^,e- rA ^[e> i+1 |^]} , (7) 

where -EfC^+ilSy is the expectation value of O at time i + 1 under the hypoth- 
esis of having the price Si at the time ti, Oj is the option value in the case 
of anticipated exercise and r is the risk-free interest. To keep under control 
the computational complexity, it is mandatory to limit the number of points 
for the integral evaluation of E[O i+1 \Si\. To this end, we can create a grid of 
possible prices, according to the dynamics of the stochastic process. Starting 
from zq, we evaluate the expectation value E'fOilSo] with p = 2m + 1, m G N 
values of z\ centered on the mean value E[zi] = z$ + A At and differing of a 
quantity of the order of a\^At 

z{ = z + AAt + jay/At, j = —m, . . . , +m. 

We can then evaluate each expectation value E[G2\z{] obtained from each one 
of the zis created above with p values for z 2 centered around the mean value 

E[z 2 \z{\ = z{ + AAt = z + 2 A At + ja VAl. 



Iterating the procedure until the maturity, we create a deterministic grid of 
points such that, at a given time ti, there are (p — l)i + 1 values of Z; L , in 
agreement with a request of linear growth. 

In the particular case of an American option, the possibility of exercise at 
any time up to the expiration date allows to develop a specific semi-analytical 
procedure, which is precise and very fast [3]. In the limit At — > and a — > 
0, the transition probability of eq. (1) approximates a delta function. This 
means that, apart from volatility effects, the price Zi at time t; L will have a 
value remarkably close to the expected value z = Zi-i + AAt, given by the 
drift. This suggests to evaluate integrals as in eq. (3) by performing a Taylor 
expansion of the kernel function Oi(e z ) around the expected value z, to arrive, 
after inserting the analytical expression (1), at the following semi-analytical 
approximation 

2 

E[Oi\Si^\ = Oi(z) + -AtO'l{z) + ..., (8) 
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Table 1 

Price of an American put option in the BSM model for the parameters t = year, 
T = 0.5 year, r = 0.1, a = 0.4, X = 10, as a function different stock prices Sq. 
The path integral 1 is performed with 200 time slices and p = 13 integration points, 
while path integral 2 is obtained with 300 time slices and p = 3. 
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where the second derivative O" can be estimated numerically. It is worth 
noticing that each expectation value E[Oi\Si_i] can be now computed once 
only p = 3 points at each slice are known, in order to evaluate the second 
derivative, provided At is taken sufficiently small. 



2.3 Numerical results and comparisons 

To test the algorithm, we present results for the particular case of the price 
of an American option in the BSM model in Tab. 1, as obtained with the 
calculation of the transition probability and the grid technique described above 
(path integral 1) and with the semi-analytical approximation of eq. (8) (path 
integral 2). As can be seen from Tab. 1, there is generally a good agreement of 
our path integral results with those known in the literature [6] and obtained by 
means of binomial trees, the finite difference method and the Green Function 
Deterministic Numerical Method (GFDNM) [6]. It is worth noticing that our 
results in the path integral 1 require only a few seconds on a Pentiumlll 
500MhZ PC, while the CPU time is negligible for path integral 2. Further 
numerical results for option prices and greek letters can be found in [3]. 



3 Neural networks 

Because the pricing formulae are given by non-linear functions and neural 
networks are known to be well suited to approximate non-linear relations, it 
is conceivable to use the fast path integral algorithm above summarized to 
generate grids of option prices over which to train learning networks. The goal 
is to construct neural networks able to price and hedge derivatives with a 
sufficient degree of accuracy to be of practical use. 



5 



Analytical 



Neural Network 




Maturity Spot Price Maturity Spot Price 

Fig. 1. Comparison for an European call option between the analytical results (left) 
and the neural network predictions (right), as a function of spot price and maturity. 
The BSM parameters are X = 10, a = 0.4, r = 0.1. 20 x 20 data points are used 
for training and 40 x 40 for generalization. 

Using the Neural Networks Toolbox of MatLab, we developed Radial Basis 
Functions networks, since it is known that this kind of networks have the best 
approximation properties for real- valued continuous functions [7]. For such 
networks, the activation functions are gaussians of variable width (a parameter 
called spread). We trained the networks using as training data the prices for 
European and American options obtained with the path integral algorithm. 
We generated several training grids of different dimensionality, starting from 
grids depending on one and two parameters to arrive at multi-parametric grids. 
After the training phase, we tested the generalization performances of the 
networks, comparing the predictions of the networks for input parameters not 
present in the training sample with the results of the benchmark. In general, 
we observed that the performances of the networks strongly depend on the 
value of the spread parameter and the number of data used in the training 
samples, as expected. 

An example of the results obtained is shown in Fig. 1, for a two-parameter 
case of an European call option in the BSM model, as a function of spot price 
and maturity. The exact analytical values (left) differ on average from the 
neural networks predictions (right) well below the 1% level. More details and 
further numerical results will be presented elsewhere [4]. 



4 Conclusions 

Computational algorithms based on the path integral approach to stochastic 
processes and neural networks can be successfully applied to the problem of 
option pricing in financial analysis. 
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The path integral can provide fast and accurate predictions for a large class of 
financial derivatives with path-dependent and early exercise features, by means 
of a careful evaluation of the transition probability and a suitable choice of 
the integration points needed to evaluate the quantities of financial interest. 
The computational cost of the algorithm is competitive with the most efficient 
strategies used in finance. Also neural networks are a powerful and flexible tool 
for option pricing. With them it is possible to parameterize functions emerging 
from financial models in a very efficient and flexible way. After an appropriate 
training, neural networks can predict option prices with good accuracy, even 
in the case of multi-parametric dependence. 

The natural development of the present work concerns the application of our 
computational methods of option pricing to more realistic models of financial 
dynamics, beyond the gaussian approximation [8]. 
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